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ABSTRACT 

We describe the generation of single-band point source catalogues from submillimetre 
HerschelSVIREi observations taken as part of the Science Demonstration Phase of 
the Herschel Multi-tiered Extragalactic Survey (HerMES). Flux densities are found 
by means of peak-finding and the fitting of a Gaussian point-response function. With 
highly-confused images, careful checks must be made on the completeness and flux 
density accuracy of the detected sources. This is done by injecting artificial sources 
into the images and analysing the resulting catalogues. Measured flux densities at 
which 50 per cent of injected sources result in good detections at (250, 350, 500) /im 
range from (fl.6, 13.2, f3.f)niJy to (25.7, 27.1, 35.8) mJy, depending on the depth 
of the observation (where a 'good' detection is taken to be one with positional offset 
less than one full-width half-maximum of the point-response function, and with the 
measured flux density within a factor of 2 of the flux density of the injected source). 
This paper acts as a reference for the 2010 July HerMES public data release. 

Key words: catalogues - submillimetre: galaxies - methods: data analysis - galaxies: 
photometry 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA, 
t E-mail: A.J.Smith@Sussex.ac.uk 



1 INTRODUCTION 

Since the discovery of the far-infrared background (FIRB; 
Puget et al. 1996; Fixsen et al. 1998; Dwek et al. 1998), suc- 
cessive surveys have aimed to identify the discrete sources 
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(primarily galaxies) responsible for this emission. With the 
launch of the ESA Herschel Space Observatory (Pilbratt 
et al. 2010), with its large (3.5 m) telescope and high sensi- 
tivity, it is now possible to resolve a much greater fraction 
of the FIRB. An essential clement of this is to have methods 
for identifying individual sources from Herschel data. 

This paper describes the generation of single-band point 
source catalogues from scan-map observations at 250, 350 
and 500 lira made using the photometer array of the SPIRE 
instrument on Herschel. The SPIRE instrument, its in-orbit 
performance, and its scientific capabilities are described by 
Griffin et al. (2010), and the SPIRE astronomical calibra- 
tion methods and accuracy are outlined by Swinyard ct al. 
(2010). The observations described here have been taken as 
part of the Herschel Multi-tiered Extragalactic Survey (Her- 
MES; Oliver et al., in preparation),^ using data from the 
Science Demonstration Phase (SDP) of the survey. These 
observations cover approximately 20 dcg^ in five regions lo- 
cated in four extra-Galactic fields, chosen for their mini- 
mal Galactic emission at far-infrared wavelengths, and for 
the amount of high-quality multi-wavelength ancillary data 
available in those fields (Oliver et al., in preparation). 

Details of the observations are given in Table 1. The 
observations in the Spitzer First Look Survey (FLS) field 
were taken in SPIRE-PACS parallel mode, at scan speed 20 
arcsec per second, while the other observations were taken 
in SPIRE-only mode, at scan speed 30 arcsec per second 
(Abcll 2218, GOODS-North and Lockman-North) or 60 arc- 
sec per second (Lockman-SWIRE). Standard SPIRE observ- 
ing modes were used for all observations. The number of 
repetitions is indicated in Table 1; for each SPIRE-only rep- 
etition, the field is scanned in both the nominal and or- 
thogonal directions, while for SPIRE-PACS parallel mode 
(FLS), one of the repetitions is in the nominal direction and 
the other is in the orthogonal direction. For the Lockman- 
SWIRE field, two separate observations were taken, offset 
from one another, in order to produce a more uniform cov- 
erage. The Abell 2218 (A2218) data were obtained through 
two observations, each consisting of 50 repetitions, separated 
by 38 days, giving complementary scan directions. All ob- 
servations were taken with nominal bias mode. More details 
are given by Oliver et al. (2010; in preparation). 

Subsets of some of the catalogues described here have 
been released to the public, as described in Appendix A. 

This paper is organised as follows. Section 2 describes 
the generation of maps and catalogues from the raw SPIRE 
data. Essential for any statistical analysis of a source cata- 
logue is to quantify the completeness and reliability of the 
catalogue, and any systematic errors in the flux densities and 
positions. This is investigated in Section 3, which describes 
a formalism for measuring these quantities and then applies 
that formalism to the catalogues. Conclusions are presented 
in Section 4. 



2 CATALOGUE GENERATION 

The data processing occurs in several distinct stages, each 
of which is described here. 



2.1 Timelines 

The SPIRE photometer contains three bolometer arrays ob- 
serving simultaneously at 250, 350 and 500 /xm. The observa- 
tions were taken as scan maps, with the telescope scanning 
the survey region at a constant rate, and with the voltage 
across each of the bolometers in the three SPIRE arrays 
being sampled at least 10 times per second (specifically, 
18.6 Hz for SPIRE-only observations and 10 Hz for obser- 
vations taken in SPIRE-PACS parallel mode, Griffin et al. 
2010). For each scan leg, this results in a series of samples 
for each bolometer, known as a 'timeline'. 

The raw timelines were processed using the standard 
SPIRE photometer pipeline (Dowell et al. 2010) to produce 
calibrated and corrected timelines in units of Jy. Specifically, 
the pipeline used was that provided in HIPE (Ott 2010) de- 
velopment version 2.0.905, with a fix applied to correct for 
a gradual drift in the astrometry (included in more recent 
versions of the pipeline), and using the following calibra- 
tion products: beam-steering mirror calibration version 2, 
fiux conversion version 2.3 and temperature drift correction 
version 2.3.2. 

A small number of cosmic ray hits ('glitches') were not 
detected by the pipeline and were propagated through to 
the maps; see Section 3.1 for a discussion of the effects this 
has on the final catalogues. 

A multiplicative correction was applied to the pipeline 
flux densities of (1.0, 1.02, 0.92) for (250, 350, 500) /xm. 
These factors were the best estimate of the correction fac- 
tors at the time the data were processed; subsequent analysis 
measured the correction factors to be (1.02, 1.05, 0.94), as 
given by Griffin et al. (2010). The current photometric ac- 
curacy of SPIRE, based on Ceres observations and models, 
is estimated to be 15 per cent (Swinyard et al. 2010) at each 
band, with a high correlation between bands. 



2.2 Maps 

From the timelines, maps were created using the default 
HIPE naive map-maker, with the default pixel sizes of (6, 10, 
14) arcsec for (250, 350, 500) /im. In map pixel i, the signal, 
di, is estimated from the Ni bolometer samples {sj} lying 
within that pixel as 



di 



N, 



1 



(1) 



while (Ti, the uncertainty in the value of di, is the standard 
error of the mean for {sj}: 



1 ' 
Ni{Ni - 1) " 



1/2 



(2) 



^ http://hermes.sussex.ac.uk 



Prior to map-making, the residual drift present in the time- 
lines, which is a residual from the temperature drift correc- 
tion (Griffin ct al. 2010), was removed by fitting a constant 
plus a linear slope to each scan timeline. Another (small) off- 
set was then applied to give the maps a mean value of zero, 
since the true (physical) zero-point for the maps is unknown. 

The overall astrometry of the maps has been adjusted 
by comparison with known positions of radio sources. This 
has typically been a correction of around 2 or 3 arcsec, which 
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Table 1. HerMES SDP SPIRE observations. For each field we give parameters for a rectangular region that avoids the edges of the 
fields, the total coverage of the observation being slightly larger. The roll angle is measured East of North for the shorter axis. (-A/samp) 
is the mean number of bolometer samples per pixel in the same typical-coverage region of the 250 /im map (6x6 arcsec pixels). The 
number of repetitions is indicated, as described in the main text. Those fields for which a Wiener filter was applied to the map data are 
indicated (see Section 2.2). 



Name 


A^rep 


RA /° 


Dec /° 


Roll /° 




Size 


samp) 


Wiener filter 


AbcU 2218 (A2218) 


100 


248.98 


66.22 


37 




9' X 9' 


1622 


N 


Spitzcr First Look Survey (FLS) 


2 


258.97 


59.39 


5 


155' 


X 135' 


30 


Y 


GOODS-North 


30 


189.23 


62.24 


42 


30 


' X 30' 


501 


N 


Lockman- North 


7 


161.50 


59.02 


1 


35 


' X 35' 


117 


N 


Lockman-SWIRE 


2 


162.00 


58.11 


2 


218' 


X 218' 


16 


Y 



is consistent with the absolute pointing error of Herschel 
(Pilbratt et al. 2010). 

Maps have been created using all of the data for each 
field, and also using two halves of the data, separated in 
time, to create two independent maps of each field, useful 
for confirmation and reliability purposes. (Note that for FLS 
one of these maps has data taken in the nominal direction 
and the other has data taken in the orthogonal direction; for 
all other fields both maps contain cross-scan data.) 

For the shallowest fields (Lockman-SWIRE and FLS), 
there was diffuse cirrus clearly visible in parts of the maps. 
In order to accentuate the signal from point sources, and 
thus to reduce the effects of the cirrus, these maps have 
been modified using a Wiener filter (Wiener 1949; Wall & 
Jenkins 2003), which is given by 



\s{fW 



\s{fr + \Nifw 



(3) 



where / is the frequency, S is the signal spectrum and A'^ 
the noise spectrum. The model for the signal is obtained 
from a noiseless simulation of sources with BLAST rmmber 
counts (Patanchon et al. 2009), which are in good agreement 
with the number counts estimated from these data (Oliver 
et al. 2010). The model for the noise is obtained from the 
difference map of the two independent maps of each field 
(which gives approximately white noise). The absolute cali- 
bration of the filtered maps is not determined at this stage; 
instead, the flux densities measured are adjusted by inject- 
ing artificial sources into the map before applying the filter 
(see Section 3.2). Fig. 1 shows part of the FLS 350 /im map 
before and after the Wiener filter has been applied. 

For these same fields, a small number of individual scans 
would have produced obvious artefacts in the final maps and 
were therefore removed. (This was due to a combination 
of steps in the thermistor timelines and the temperature 
drift correction used.) In these regions, the coverage has 
consequently been reduced by approximately 25 per cent. 
For FLS, out of 117 scans, this has affected two scans ob- 
tained at 250 fim and one scan obtained at 500 /J,m, while for 
Lockman-SWIRE, out of 160 scans, the number affected at 
(250, 350, 500) /im was (3, 1, 2), respectively. 



2.3 Catalogues 

Source catalogues have been generated for each band in each 
field. Details of the columns are given in Appendix A. Source 




Figure 1. Part of the 350 /^m map from FLS, covering an area 
of approximately 36 arcmin by 33 arcmin, before (left) and af- 
ter (right) applying the Wiener filter. The point sources are sig- 
nificantly emphasised by the Wiener filter, thus diminishing the 
efi'ects of the large-scale diffuse cirrus. 



flux densities have been estimated using the SUSSEXTrac- 
TOR point source extractor (Savage & Oliver 2007) as im- 
plemented in HIPE 3.0. For computational efficiency, a Gaus- 
sian point-response function (PRF) was assumed, with full- 
width half-maximum (FWHM) of (18.15, 25.15, 36.3) arcsec 
for (250, 350, 500) /xm, and with a Gaussian approximation 
used for the beam area: 



7r(FWHM)^ 
4 In 2 



(4) 



Griffin et al. (2010) state that Gaussian beams with FWHM 

of (18.1, 25.2, 36.6) arcsec provide a good approximation to 
the true beam, so we assume the errors introduced by our 
choice of PRF will be small compared with other sources of 
uncertainty. 

The ffux density is given by 



djVi 

^2 



V? 



(5) 



where the summation is over a local region around the source 

position, di is the value of the map pixel, ai is the value of 
the error map pixel and Vi is the smoothing kernel (matched 
filter). The uncertainty in the ffux density is discussed in 
Section 2.4. 

For isolated point sources and white noise, the optimal 
matched filter is the beam itself. But for a higher density of 
sources the optimal filter will be narrower than the beam, 
up to the limit of complete confusion (no white noise), in 
which the optimal filter is a deconvolution filter, given by 
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the inverse of the beam (see Chapin ct al. 2010, appendix 
A). 

However, the method used here, for convenience, was to 
take as the smoothing kernel the central region of a Gaus- 
sian PRF centred on the pixel closest to the source position. 
SUSSEXTRACTOR was used in two distinct ways, one for the 
shallower fields, and one for the deeper, so as to deal with 
the high pixel-to-pixcl noise in the shallower maps and to 
exploit the high signal-to-noise in the deeper maps. 

For the shallower fields (FLS, Lockman-SWIRE and 
Lockman-North) , a single pass was performed, which in- 
volves SUSSEXTRACTOR applying a smoothing kernel to the 
image (the central 5x5 pixels of the Gaussian PRF), to 
obtain an image in which the value of each pixel is the max- 
imum likelihood estimate of the flux density of a source cen- 
tred on the centre of that pixel, as given in Equation (5). 
In this smoothed map, the extractor then searched for lo- 
cal maxima, comparing each pixel with its eight immediate 
neighbours, and the value at these peak positions was taken 
as the estimate of the source flux density. The position of 
the source was refined to sub-pixel accuracy based on the 
intensity of the surrounding pixels. 

For the deeper fields (A2218 and GOODS- North), first 
SUSSEXTRACTOR was run using no smoothing, in order to 
find the positions of the local maxima in the image. Then 
SUSSEXTRACTOR was run again, this time with a small 
smoothing kernel (the central 3x3 pixels of the Gaussian 
PRF), in order to estimate the source flux densities at these 
position. This method has been adopted in order to extract 
to fainter flux densities, and in order to reduce the number 
of close pairs of sources that are blended into a single source 
in the catalogue. 

Sources have been extracted close to the edge of the 
images. However, a central region has been defined for each 
field, so that sources can be selected within a simple rect- 
angular region of the image, avoiding the edges. Parameters 
defining these regions arc in Table 1. These regions have 
an easily-determined area, and therefore the subsets can be 
used for studies of the number density of sources. 

At the positions of the sources in the catalogues, flux 
densities have been estimated from the two independent 
maps, each produced from half of the data (see Section 
2.2). These flux densities have been included in the cata- 
logues, and may be used for investigations of the reliability 
of sources (see Section 3.1). 

A multiplicative factor has been applied to all flux den- 
sities, in order to give approximately zero mean offset in 
log-fiux density for the brightest injected sources (see Sec- 
tion 3.3). This is to account for the arbitrary normalisa- 
tion and other effects of the Wiener filter, and also for a 
systematic underestimation of the flux densities by SUSSEX- 
TRACTOR, primarily due to the assumption in the flux den- 
sity estimation that the source centre is aligned with the 
centre of a pixel. First, some bright sources of flux den- 
sity S were added to the images. Then the flux densities 
of these sources were measured. If the mean measured flux 
density was S'mean, then the multiplicative factor was cho- 
sen to be S/Smean- FoT the (250, 350, 500) fim flux densities, 
the factors applied to the standard (naive map) data were 
A2218: (1.052, 1.062, 1.040), GOODS-North: (1.067, 1.062, 
1.074), Lockman-North: (1.028, 1.040, 1.038). For FLS and 
Lockman-SWIRE, the factors applied to the Wiener-filtered 



data were (1.266, 1.304, 1.482) for FLS and (1.613, 1.580, 
1.713) for Lockman-SWIRE. (Note that these factors were 
derived before certain improvements were made to the 
method of injecting artificial sources, so the offset in log- 
flux density is only approximately zero for bright injected 



2.4 Uncertainties in source flux densities 

The formal uncertainty in the fiux density in Equation (5) 
is given by 



as = 1 



(6) 



It should be noted that in expressing the uncertainty in the 
flux in this way, it has been assumed implicitly that the co- 
variance between pixels is negligible. For the naive maps, 
this is a reasonable assumption, but for the Wiener-filtered 
data (FLS and Lockman-SWIRE) the covariance between 
neighbouring pixels is not negligible, and thus the formal 
uncertainty in the flux density from Equation (6) will be a 
poor estimate of the formal uncertainty. Moreover, for the 
Wiener-filtered data, the error map (which was not Wiener 
filtered) has also been scaled by the (large) multiplicative 
factors given in Section 2.3, which will significantly increase 
the value of as- These values should therefore not be over- 
interpreted for the Wiener-filtered data. However, the to- 
tal noise estimates (see below) are based primarily on the 
statistics of the smoothed images provided by the source ex- 
traction software, rather than on the formal uncertainty in 
Equation (6), so these are more robust against these effects. 

These flux density uncertainties are believed (for the 
non- Wiener-filtered data) to give a fair estimate of the in- 
strumental noise, and will be referred to as such hereafter. 
But they do not include the effects of source confusion (the 
high density of sources relative to the size of the SPIRE 
beams) nor the effects of the uncertainty in the PRF or the 
SPIRE absolute flux calibration. So the true uncertainty in 
the source flux density will be significantly higher than the 
instrumental noise. 

The total noise, taking account of confusion noise as 
well as instrumental noise, is estimated as follows: 

(i) The smoothed map is obtained from the source ex- 
traction software, with the value in each pixel, j, being an 
estimate of the flux density, Sj, of a point source assumed 
to lie at the centre of that pixel, calculated from Equation 
(5). 

(ii) The typical total noise in the source flux density, 
ftotai, is derived from the statistics of this smoothed map: 



""total — 



{Sj - median(g,))^ 
N 



(7) 



where the summation is over the central region of the map, 
restricted to those pixels in which Sj < median(5,), and 
where N is the number of pixels included in the summation. 

(iii) The typical instrumental noise, iTinstrumentai, is esti- 
mated as the median of Equation (6), over the central region 
of the map: 



f instrumental = median (crs,; 



(8) 
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Table 2. Approximated la uncertainty in the flux density of a 
typical point source, in mjy, from the combined effects of in- 
strumental and confusion noise, as described in the text. Shown 
in parentheses is the median Icr instrumental noise in the flux 

density measurement of a point source, in mJy. For the Wiener- 
filtered data (FLS and Lockman-SWIRE), the instrumental noise 
is believed to be over-estimated, and is shown in italics. 



Field 



f total (^'instrumental)! rn^Y 

250 /im 350 ixm 500 ixva 



A2218 


5.9 


(0.6) 


7.6 


(0.6) 


7.8 


(0.6) 


FLS2 


8.8 


(3.1) 


10.0 


(3.2) 


11.1 


(4.5) 


GOODS-North 


5.7 


(0.9) 


7.4 


(0.9) 


7.8 


(1.1) 


Lock.-North 


7.0 


(1.7) 


8.5 


(1.7) 


8.8 


(2.1) 


Lock.-SWIRE2 


10.4 


(6.6) 


11.6 


(6.5) 


11.8 


(8.7) 



(iv) A single value for the typical contribution from 
source confusion to the uncertainty in the flux density of 
a point source is derived by subtracting, in quadrature, the 
typical instrumental noise from the typical total noise: 



^confusion 



2 

"■total 



"instrumental • 



(9) 



This is a similar, but not identical, quantity to the 'con- 
fusion noise' derived by Nguyen et al. (2010), which is the 
contribution from source confusion to the uncertainty in the 
intensity of a typical map pixel and is measured to be '5.8, 
6.3 and 6.8 mjy/beam at 250, 350 and 500 //m, respectively'. 

(v) The total noise for each individual source (pixel j) is 
then obtained by adding, in quadrature, this typical value for 
the confusion noise to the source's own instrumental noise 
value, given by Equation (6): 

"instrumental, J "confusion ^^ "S.J ■ (-^^) 

The values for the typical total noise and the typical in- 
strumental noise are shown in Table 2. Note that smoothing 
of the maps increases the confusion noise (although it de- 
creases the instrumental noise contribution), since smooth- 
ing increases the size of the effective beam. More smoothing 
has been applied to the shallower fields (FLS, Lockman- 
North, Lockman-SWIRE) than to the deeper fields (A2218, 
GOODS-North), so greater confusion noise is to be expected 
in those fields. Moreover, the FLS and Lockman-SWIRE 
fields have been smoothed using a Wiener filter (with differ- 
ent filters for each field), which will broaden the PRE and 
thus increase the confusion further. 

As discussed above, the measurement of the instrumen- 
tal noise for the Wiener-filtered data (ELS and Lockman- 
SWIRE) is believed to be over-estimated, so it is not possible 
to give a reliable estimate of the confusion noise using Equa- 
tion (9). However, for the remaining data, (Tconfusion is found 
to be (5.9, 7.5, 7.7) mJy for A2218, (5.6, 7.4, 7.7) mJy for 
GOODS-North and (6.8, 8.3, 8.5) mJy for Lockman-North, 
all for (250, 350, 500) /im respectively. 

The initial threshold on the catalogues is 3(t, based on 
the instrumental noise. Some further cuts have been applied 
to the released catalogues (see Appendix A). 



^ Estimates of the noise for FLS and Lockman-SWIRE were re- 
vised slightly since the publication of Schulz et al. (2010), but the 
findings of that paper are unaffected. 



The raw source counts for the central region of each 
field are shown in Eig. 2 with the total number of sources in 
the central region of each field, and the number with signal- 
to- (total) noise ratio greater than 3, shown in Table 3. 



3 QUALITY OF THE CATALOGUES 

A catalogue will be of limited use without some measure or 

assurance of its quality. This may be the reliability (number 
of false detections) , completeness (probability that a genuine 
source will be included in the catalogue) or the accuracy of 
the parameters of the sources (position and fiux density). 
These will be discussed below. 



3.1 Reliability 

The reliability of a source catalogue is conventionally a mea- 
sure of the fraction of detections, at a given flux density, that 
are spurious. A spurious detection may happen as a result of 
noise in the map pixels (due to a small number of bolometer 
samples, each with a significant uncertainty), or as a result 
of other factors contributing to the detector signal, such as 
any cosmic ray hits (glitches) that are not removed by the 
pipeline. 

When the noise in the data is due entirely to these (in- 
strumental) effects, the probability that a detection is gen- 
uine (or spurious) can be estimated from the signal-to-noise 
ratio of the source. However, in these _ffersc/ie/-SPIRE data, 
the dominant source of noise is confusion, that is, the mea- 
surement of the fiux density of any particular source being 
contaminated by the flux density of neighbouring sources. 
This means that the signal-to- (total) noise of a detection 
cannot be used in any straightforward way to give the prob- 
ability that it is spurious. 

The number of such spurious detections that would arise 
from instrumental noise may be estimated using the maps 
and catalogues generated from the two halves of the data 
for each field (see Section 2.2). 

If di is the measured intensity in a map pixel from the 
first half of the data and 6,2 is the measured intensity in the 
same pixel from the second half, then the intensity for the 
total map will be 

dr + d, _ ^^^^ 



'total — 

Two 'difference maps' may be obtained by taking the differ- 
ence between these two measurements: 

d2' 



t^diiference ■ 



(12) 



This difference map is then an instrumental-noise map with 
astronomical fiux (and confusion noise) removed. Executing 
the source extraction on these maps will give an estimate of 
the number of spurious detections that might be expected, 
in the absence of confusion noise. 

Any unremoved cosmic ray hits will either leave a pos- 
itive spike or a negative spike in the difference map, de- 
pending on which half of the data is affected. The source 
extraction is therefore executed on both the 'positive' and 
'negative' difference maps, from Equation (12). 

Table 4 shows the numbers of sources detected from 
these difference maps with additional details given below. 
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Flux density, mjy 

[ O O SOO^m ■ 350/; m C O 250y<m | 

(e) Lockman-SWIRE 

Figure 2. Raw source counts for the five fields, showing the number of sources per unit area per logj^Q interval in the central region. No 
corrections have been applied for incompleteness or fiux boosting. Error bars are the counts in that bin divided by the square root of the 
number of sources in the bin. Each plot shows the results for 250 /.tm (blue), 350 /^m (green) and 500 /.im (red). 
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Table 3. Number of sources in the central region of each field with signal-to- (total) noise greater than 3. In parentheses are the number 
of sources with signal-to-(instrumental) noise greater than 3. 



Field Number of sources > 3crtotal (> 3crinstrumental) 

250 /im 350 /im 500 fira 



A2218 


41 


(119) 


12 


(64) 


5 


(36) 


FLS 


3946 


(12 862) 


1822 


(7120) 


637 


(2751) 


GOODS-North 


385 


(1421) 


150 


(713) 


48 


(344) 


Lockman-North 


325 


(1082) 


141 


(586) 


61 


(255) 


Lockman-SWIRE 


6731 


(13 890) 


2757 


(7867) 


836 


(1902) 



Table 4. Number of detections in the central regions of the difference maps greater than 3(Ttotal- Each pair of numbers is the number of 
detections in the first and then the second difference map, corresponding to the positive and negative forms of Equation (12), respectively. 
In parentheses are the number of detections with flux density greater than 3<Tinstrumental • The numbers should be compared with those 
given in Table 3. 



Field Sources > 3(Ttotal (> 3CTinstrumental) 

250 Mm 350 /im 500 //m 



A2218 


0-1-0 


(93-1-108) 


0+0 


(62+66) 


0+0 


(35+23) 


FLS 


7-1-3 


(28-1-25) 


4+4 


(11+10) 


0+1 


(4+6) 


GOODS-North 


04-0 


(2+1) 


0+0 


(3+0) 


0+0 


(1+0) 


Lockman-North 


0+0 


(9+6) 


0+0 


(5+7) 


0+0 


(5+36) 


Lockman-SWIRE 


2-1-6 


(3+7) 


1+1 


(1+1) 


4+3 


(4+4) 



The numbers of detections should bo compared with those 
in Table 3. 

For A2218 (0.022 deg^), all of the 'detections' in the 
difference maps have flux densities below 7.2 mjy. The de- 
tections arc generally found along stripes in the map that re- 
main as a result of our method for subtracting baselines from 
the timelines (Section 2.1). A small number of detections in 
the difference maps are associated with bright sources: these 
may have arisen as a result of the strong gradients in the sig- 
nal associated with the sides of the beam, or as a result of 
the ellipticity of the beam (7-12 per cent, Griffin et al. 2010) 
and the change in position angle between the two observa- 
tions. 

For FLS (5.8deg2), 7+25 of the 'detections' at 250 /im, 
5+9 at 350 /im and 0+6 at 500 /im are above the la total 
noise values in Table 2, with one of these having a flux den- 
sity greater than 200mJy. For Lockman-SWIRE (13deg^), 
the measured flux densities arc between 28 and 400 mJy. By 
inspection of the images, all of these detections were found 
to be due to glitches that were not removed by the pipeline, 
with the exception of 12 detections in FLS at 250 /im: 4 as- 
sociated with bright sources, and 8 detections fainter than 
10.6 mJy. 

For GOODS-North (0.25 deg^) and for Lockman-North 
(0.34deg^), all of the measured flux densities are below 
7.5 mJy, with the exception of the second difference map for 
Lockman-North at 500 /tm, which gives 33 detections along 
one scan line, due to residuals in the baseline subtraction, 
having flux densities between 6.3 and 12.7 mJy. 

The spurious detections due to unremoved glitches 
cause a detection in the map from one half of the data but 
not in the map from the other half. The flux densities gen- 
erated from these half-data maps are included in the cata- 
logues, and may be used to identify some of these. 



3.2 Completeness, flux density accuracy and 
positional accuracy: method 

The completeness, flux density accuracy and positional ac- 
curacy have been investigated by injecting artificial sources 
into the timelines, and then creating new maps from those 
modified timcfincs. This ensures that both the signal and its 
uncertainty are modified by the injected sources, as given by 
Equations (1) and (2). For the Wiener-filtered data, the er- 
ror maps were created in this way, but the signal maps were 
created by injecting sources directly into the unfiltered im- 
ages, and then applying the Wiener filter to the images. 

For each field, this is done multiple times, with the same 
flux density for all of the injected sources. Sources arc placed 
on a grid, with spacing offset from the pixel size of the im- 
age and large enough that the sources can be treated in- 
dependently (as an approximation to the process of adding 
sources one at at time at random positions). (The precise 
spacing chosen was 113.387 arcsec for A2218, GOODS-North 
and Lockman-North, increased by a factor of 3 for FLS and 
Lockman-SWIRE.) An example is shown in Fig. 3. 

The procedure for measuring the completeness and flux 
density accuracy is as follows. 

First the source catalogues are produced: 

(i) The source extraction is performed on the map with 
no artiflcial sources added, as described in Section 2.3, in 
order to define a reference catalogue for each band. 

(ii) A truth catalogue is created, consisting of the grid 
of artificial sources, and this is used to create maps with 
injected sources, with all such sources having the same flux 
density. The whole procedure is repeated with each iteration 
having a different flux density for the injected sources. The 
flux densities chosen are 1, 2, 3, 4, 7, 10, 20, 30, 40, 70, 100, 
200, 300, 400, 700, 1000 and 4000 mJy. 



8 A.J. Smith et al. 




Figure 3. The 350 /im map from GOODS-North, with a grid of 
100 mJy sources injected. 



(ill) For each map with injected sources, the source ex- 
traction is performed, in exactly the same way as for the ref- 
erence catalogue, to produce additional source catalogues, 
to be compared with the reference catalogues. 

Next, with the reference catalogue, truth catalogue of 
artificial sources, and source catalogues for each injected flux 
density, the catalogues are compared as follows, for each 
band and for each injected flux density: 

(i) The artiflcial source truth catalogue is flrst cross- 
matched with the reference catalogue from the real data. 
For each source in the truth catalogue, the closest match 
within 1 times the FWHM is chosen (if such a source is 
present). If this match has a flux density within a factor of 
2 of the injected source flux density, the match is identified 
as a 'good' match. Any such 'good' matches are discarded 
from further analysis; otherwise, when these (serendipitous) 
matches are included, the measured completeness can be 
misleading, particularly for source extraction methods that 
produce a large number of spurious, faint detections. How- 
ever, excluding these sources will have a small effect on the 
estimates of the completeness, because part of the incom- 
pleteness comes from the fact that sources can be too close 
to other sources and therefore not counted. 

(ii) Next the truth catalogue (without the serendipitous 
sources from the previous step) is cross-matched with the 
source catalogue derived from the map with injected sources. 
'Good' matches are found, as above. The completeness is de- 
fined as the number of good matches divided by the number 
of injected sources (minus the serendipitous sources). 

For example, if 200 sources with fiux density 30 mJy are 
injected into the map, but 20 of those already (by chance) 
have 'good' counterparts in the original map (without in- 
jected sources), then the remaining number of sources is 180. 
If 162 of these have good matches in the source list extracted 
from the map with injected sources, then the completeness 
at 30 mJy is 162/180 = 90 per cent. 

(iii) The fiux density and positional accuracy are found 
by comparing the extracted fiux densities and positions with 
the injected fiux densities and positions. 



Table 5. Flux density corresponding to 50 per cent completeness 
in Fig. 4. 



Field 


Flux density. 


mJy 




250 Mm 


350 /^m 


500 /im 


A2218 


13.4 


15.9 


15.5 


FLS 


17.4 


19.0 


22.1 


GOODS-North 


11.6 


13.2 


13.1 


Lockman-North 


13.6 


15.7 


17.0 


Lockman-SWIRE 


25.7 


27.1 


35.8 



3.3 Completeness, flux density accuracy and 
positional accuracy: results 

The completeness, as defined above, for each field is shown in 
Fig. 4. Treating the completeness, C, as the parameter of a 
binomial distribution, the posterior probability for the value 
of the completeness being C is given by a beta distribution: 

P(Ci7Vi„j, 7V,ec) oc C^'-(l - C)^'"^-^'- (13) 

where A'^inj is the number of injected sources and A^roc is the 
number of sources recovered. This is used to obtain the error 
bars on the completeness in Fig. 4. 

For each completeness curve, the fiux density is found 
at which the completeness is 50 per cent; these values are 
given in Table 5. 

The accuracy of the fiux densities and positions of the 
recovered sources are shown in Figs. 5 and 6 respectively. 

Table 6 shows the mean offset in flux density for various 
input flux densities, while Table 7 shows the corresponding 
RMS scatter in the flux density offset and Table 8 shows the 
RMS scatter for the offset in RA. 

Several points should be noted from these. 

(i) The source extraction method has been adjusted by 
applying a multiplicative factor to all flux densities in or- 
der to give good recovered flux densities for bright injected 
sources. These factors are given in Section 2.3. This can be 
seen in Fig. 5 by the way the flux density offset is mea- 
sured to be approximately zero for bright injected sources 
for most of the fields (some late adjustments were made to 
the Wiener filtered data, leaving a very small residual offset 
for bright flux densities for FLS and Lockman-SWIRE). 

(ii) In Fig. 5, towards fainter flux densities, there is a 
trend of a bias towards an overestimated flux density and 
decreasing RMS scatter; this is a selection effect due to flux 
boosting and the requirement that a good match will have 
a flux density within a factor of 2 (0.3 dex) of the input flux 
density. 

(iii) Also in Fig. 5 there is a dip in the plots at around 
10-40 mJy, corresponding to an underestimate of the flux 
density. This is likely to be due to the maps having a zero 
mean, leading to a systematic underestimation of the flux 
densities. Corrective factors were applied to the flux densi- 
ties in order to give good agreement for bright input sources, 
but these were multiplicative corrections, rather than correc- 
tions with both an additive and multiplicative component, 
and thus the corrections have been effective only for bright 
input sources, where the multiplicative factor is dominant. 

(iv) The scatter in the recovered flux densities in Fig. 5 
may be compared with the uncertainties in the flux densities 
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Input flux density (mJy) Input flux density (mJy) 



SOO^m : ^ 350/fni Q Q 23n,«m Q Q .TOO/m " 350/fni Q Q l^Qftm 

(c) GOODS-North (d) Lockman-North 




Input flux density (mJy) 

I O O Sm^m Q O 35Qam O O 250/;m | 

(e) Lockman-SWIRE 

Figure 4. Completeness, for the five fields, as in Fig. 2. The flux density corresponding to 50 per cent completeness is given in Table 5. 
Error bars and dashed lines are the lower and upper bounds of the posterior probability distribution for the completeness, chosen such 
that the probabilities at the bounds are equal, and so that the integrated probability between the bounds is 68 per cent. 
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Figure 5. Flux density accuracy, for the five fields, as in Fig. 2. Error bars are the RMS of the logio (recovered flux density / input flux 
density) at that flux density. 
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Figure 6. Positional accuracy in RA (offset to the east, in arcsec), for the five fields, as in Fig. 2. Error bars are the RMS of the positional 
error at that flux density. Similar results are found for the offset in declination. 
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Table 6. Mean value of logio (recovered flux density / Input flux density), for each band, at Input flux densities of 20, 40, 100 and 
1000 mjy, from Fig. 5. 



Field Mean offset, dex 

250 iim 350 /im 500 /im 





20 


40 


100 


1000 


20 


40 


100 


1000 


20 


40 


100 


1000 mJy 


A2218 


-0.086 


-0.013 


-0.001 


0.013 


-0.035 


-0.018 


0.008 


-0.003 


-0.042 


-0.010 


0.005 


0.003 


FLS 


-0.040 


-0.049 


-0.028 


-0.039 


-0.021 


-0.044 


-0.022 


-0.038 


0.074 


-0.046 


-0.037 


-0.049 


GOODSN 


-0.057 


-0.008 


0.004 


0.013 


-0.038 


-0.006 


0.004 


-0.003 


-0.010 


0.012 


0.015 


0.006 


LOCKN 


-0.044 


-0.010 


0.002 


-0.006 


-0.021 


-0.008 


0.008 


-0.015 


0.021 


0.015 


0.016 


-0.008 


LOCKSW 


0.132 


-0.034 


-0.012 


-0.020 


0.126 


-0.054 


-0.051 


-0.050 


0.225 


0.007 


-0.056 


-0.054 



Table 7. RMS scatter in logio (recovered flux density / input flux density), for each band, at input flux densities of 20, 40, 100 and 
1000 mJy, from Fig. 5. 



Field RMS scatter, dex 







250 








350 


/im 






500 /im 






20 


40 


100 


1000 


20 


40 


100 


1000 


20 


40 


100 


1000 mJy 


A2218 


0.125 


0.066 


0.038 


0.014 


0.154 


0.085 


0.055 


0.030 


0.118 


0.077 


0.034 


0.017 


FLS 


0.134 


0.092 


0.044 


0.010 


0.139 


0.107 


0.050 


0.010 


0.103 


0.112 


0.057 


0.017 


GOODS-North 


0.110 


0.068 


0.039 


0.018 


0.127 


0.079 


0.042 


0.027 


0.126 


0.079 


0.047 


0.026 


LOCK-North 


0.117 


0.070 


0.036 


0.017 


0.136 


0.083 


0.043 


0.025 


0.141 


0.088 


0.047 


0.025 


LOCK-SWIRE 


0.082 


0.104 


0.049 


0.014 


0.083 


0.111 


0.059 


0.018 


0.047 


0.087 


0.063 


0.020 



given in the catalogues. The latter are, to first order, inde- 
pendent of the flux density of the source, being based on the 
statistics of the whole map (see Section 2.4). However, for 
bright injected sources, the former is approximately propor- 
tional to the input flux density. This is because the dominant 
source of noise in the measurement of the flux density of a 
bright source is not confusion noise but rather systematic 
errors from the source extraction method, for example, vari- 
ations in the recovered flux density depending on where the 
centre of a source lies within a pixel (sec Section 2.3). For 
bright sources, therefore, the values given in Table 7 should 
be used as an approximate guide to the uncertainty in the 
measured flux. 

(v) The quality of the positions and flux densities for the 
brightest flux densities can be seen from Tables 7 and 8 
to depend on the approach used to extract the sources. 
This is discussed in Section 2.3 and depends on whether 
the shallow approach has been used (FLS, Lockman-North 
and Lockman-SWIRE), or the deep approach (A2218 and 
GOODS-North) and whether a Wiener filter has been ap- 
plied to the maps (FLS and Lockman-SWIRE). In partic- 
ular, where a smaller amount of smoothing has been ap- 
plied to the maps in order to improve the extraction of faint 
sources and deal with the problem of source blending (A2218 
and GOODS-North), the scatter in the recovered flux den- 
sities and positions for bright injected sources is larger than 
when more smoothing is applied. 

(vi) The combined effect of the flux density uncertain- 
ties in Fig. 5 and the steep number counts seen at SPIRE 
wavelengths (Oliver et al. 2010) will mean that at any given 
measured flux density, most of the sources will have a true 
flux density which is fainter than the measured flux density, 
even if there is no systematic offset of injected to measured 
flux densities, as is the case here. This phenomenon ('flux 



boosting') must be taken into account when estimating the 
true flux densities of sources. See Oliver et al. (2010) for 
further discussion. 

(vii) Some features may be discerned for bright injected 
flux densities in Fig. 5, such as a departure away from a 
horizontal slope. This is due to the effect of the injected 
sources on the error maps. In the steep slope of the point- 
response function, there is a large scatter in the intensity 
of the bolometer samples falling within the map pixel. This 
leads through Equation (2) to a higher value for the uncer- 
tainty for that map pixel, which leads through Equation (6) 
to that map pixel being given a lower weight in the source 
extraction. These changes in the relative weight given to the 
pixels in a point source lead to changing estimates of the 
source flux density, based on the flux density of the injected 
source. 

(viii) Finally, it should be noted that the method used 

here to evaluate the accuracy of the measured flux densities 
and positions uses idealised artificial sources, which will be 
subtly different from real sources. One cause of these differ- 
ences would be the assumed Gaussian PRE; the Airy rings 
around an extremely bright source would influence the de- 
tections close to that source. For the very brightest sources, 
the best way to test SUSSEXTRACTOR (and other algorithms) 
would be on real observations of SPIRE calibration sources 
(see Swinyard et al. 2010). 



4 CONCLUSIONS 

We have described the approach adopted to generate single- 
band catalogues from HerMES SDP SPIRE observations, 
some of which have been made publicly available. A for- 
malism has been developed to assess the quality of these 
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Table 8. RMS scatter in the offset in RA, for each band, at input flux densities of 20, 40, 100 and 1000 mjy, from Fig. 6. 
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Field RMS scatter, arcsec 







250 








350 








500 /im 






20 


40 


100 


1000 


20 


40 


100 


1000 


20 


40 


100 


1000 mJy 


A2218 


2.821 


1.533 


1.019 


0.727 


6.337 


2.301 


1.466 


1.491 


4.980 


3.200 


1.631 


1.472 


FLS 


4.935 


2.418 


0.976 


0.403 


7.679 


3.954 


1.528 


0.360 


12.272 


6.628 


2.440 


0.668 


GOODS-North 


2.439 


1.564 


1.283 


1.035 


4.030 


2.546 


1.554 


1.681 


5.883 


3.418 


2.604 


2.420 


LOCK-North 


3.196 


1.746 


0.795 


0.328 


5.935 


2.846 


1.396 


0.815 


9.707 


5.206 


1.993 


0.958 


LOCK-SWIRE 


5.428 


2.971 


1.190 


0.681 


9.077 


5.552 


2.258 


1.166 


13.314 


7.916 


2.922 


1.226 



catalogues, and recommendations have been made for usage 
of the catalogues based on these results. 

Possible improvements to the method used here have 
been identified above, including the following: 

(i) The filtering of the data could be refined, both by 
including information about cirrus in the noise spectrum 
for the Wiener filter (Section 2.2), and by using optimized 
matched filters for the detection and measurement of point 
sources (Section 2.3). 

(ii) A better estimation of the background would deal 
with the additive offset to the measured flux densities, dis- 
cussed in Section 3.3. This could be achieved either by deter- 
mining a physical zero-point for the maps or through local 
background estimation for each source. 

(iii) A better estimation of the fiux densities could be 
achieved by relaxing the implicit assumption that all sources 
lie in the centre of map pixels. This could be achieved by 
first finding the source positions to sub-pixel accuracy, as 
in the current method, and then choosing an appropriate 
smoothing kernel in Equation (5) based on this measured 
position. 

Only one approach has been presented here for source 
extraction. Other approaches exist, both in terms of other 
algorithms, and alternative ways of using any particular al- 
gorithm, such as an iterative approach, removing the bright- 
est sources from the image at each iteration. Alternatively, 
information from multiple bands may be used, either simul- 
taneously, to extract sources in multiple Herschel bands at 
the same time, or by using prior positions from other wave- 
lengths, a method used for many HerMES results and de- 
scribed by Roseboom et al. (2010). Future work within Her- 
MES will explore these approaches in more depth. 
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APPENDIX A: DATA RELEASE 

The central region of the maps for A2218 and part of the 
catalogues for A2218, Lockman-SWIRE and FLS described 
in this paper have been released to the public on the Herschel 
Database in Marseille (HeDaM).^ 

For A2218, three single-band SPIRE catalogues have 
been released, including all sources within the central re- 
gion with flux densities greater than 20mJy (6 sources at 
500 /um, 19 at at 350 /xm and 35 at 250 /um). Descriptions of 
the columns are given in Table Al. 

For FLS and Lockman-SWIRE, 250 /im catalogues have 
been released, containing all sources within the central re- 
gion with flux densities greater than 100 mjy, as long as the 
source has one and only one counterpart within 10 arcsec in 
the associated 24 /xm catalogue. The resulting FLS catalogue 
contains 45 sources and the Lockman-SWIRE catalogues 
contains 114. For FLS and Lockman-SWIRE, in addition 
to the columns described in Table Al, additional columns 
are provided in the catalogues, some derived from a fusion 
of ancillary data (Vaccari et al., in preparation) and others 
containing SPIRE list-driven quantities derived from 24 fmi 
source positions Roseboom et al. (2010). 



^ http://hedam.oamp.fr/HerMES/ 



HerMES: Herschel-SPIRE point source catalogues 

Table Al. Names and descriptions of the columns included in the public data release (see Appendix A). 



Name 


Col. No. 


Description 


name 


1 


HerMtb ID 


ra 


2 


Right Ascension (deg) 


dec 


3 


Declination (dcg) 


raErr 




Right ascension uncertainty (dcg) 


decErr 


5 


Declination uncertainty (dcg) 


flux 


6 


Source flux density (mjy) 


fluxErr 


7 


Source flux density formal uncertainty (instrumental noise, mJy) 


quality 


8 


Signal to instrumental noise: flux/fluxErr 


index 


9 


Sequential number in full catalogue 


centralRegion 


10 


True if source lies within a well defined central region of the map 


fluxErrTotal 


11 


Total uncertainty in the source flux density, due to confusion and instrumental noise (mJy) 


SNR 


12 


Signal to total noise: flux/fluxErrTotal 


fluxHalfDatal 


13 


Source flux density, as measured using a map based on the first half of the data (mJy) 


qualityHalfDatal 


14 


Signal to instrumental noise, as measured using a map based on the first half of the data 


fluxHalfData2 


15 


Source flux density, as measured using a map based on the second half of the data (mJy) 


qualityHalfData2 


16 


Signal to instrumental noise, as measured using a map based on the second half of the data 


extended 


177/178 


If true, source noticeably extended in the SPIRE 250 /um image (flux density should be treated 



with caution) [not A2218] 



